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Abstract. We obtain exact periodic solutions of the positive and negative 
modified Kortweg-de Vries (mKdV) equations. We examine the dynamical 
stability of these solitary wave lattices through direct numerical simulations. 
While the positive mKdV breather lattice solutions are found to be unstable, the 
two-soliton lattice solution of the same equation is found to be stable. Similarly, 
a negative mKdV lattice solution is found to be stable. We also touch upon the 
implications of these results for the KdV equation. 
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1. Introduction 

Among the soliton bearing nonlinear integrable equations, sine-Gordon, nonlinear 
Schrodinger (NLS), Korteweg-de Vries (KdV) and the modified Korteweg-de Vries 
(mKdV) are of special interest |3 El |S] . These equations possess exact breather 
solutions. Therefore, they may also have exact solutions in the form of a spatially 
periodic array of single breathers, i.e. breather lattices. Similarly to the sine-Gordon, 
NLS and KdV, the mKdV equation is also important in many physical contexts. 
For example, it appears in the context of ion acoustic solitons jS], van Alfven waves 
in collisionless plasma T, Schottky barrier transmission lines [8], models of traffic 
congestion 9 as well as phonons in anharmonic lattices |1U| . Furthermore, the 
modelling of a subclass of hyperbolic surfaces [HI , slag-metallic bath interfaces as 
well as meandering ocean jets pi] is also related to the mKdV equation. The dynamics 
of thin elastic rods has also been demonstrated to be reducible to the mKdV equation 
|14| . Finally, if one studies the examples of surface dynamics that are purely local, 
yet maintain global constraints like conservation of perimeter and enclosed area, one 
finds that these dynamics are closely related to the KdV and mKdV hierarchies |15| . 

The nonlinear term in the mKdV equation (Qu 2 u x ) may assume either a positive 
or a negative sign. We will classify the equation as "positive mKdV" and "negative 
mKdV" if the prefactor is +1 or —1, respectively. In an earlier study, the stability of 
the sine-Gordon breather lattice was examined Recently, a particular form of an 
exact breather lattice solution was obtained in ^7] for the positive mKdV equation 
and its stability was discussed. 

The aim of the present Letter is to present a new class of periodic solutions both 
for the positive and for the negative mKdV equations. We also intend to examine 
the dynamical stability of these novel classes of solutions and particularly to illustrate 
that many of them can be dynamically stable. This is notably different from the 
behavior observed previously for breather lattice solutions in the models of |161 [T7| . 
Furthermore, it is worth noting that the solutions of the negative mKdV are related 
to those of the KdV equation via the Miura transform ^B]- Thus, our solutions can 
be translated into exact periodic solutions of the KdV equation as well. The latter is 
also a ubiquitous equation in a variety of fields ranging from conformal field theory to 
plasma physics. 

Our presentation is structured as follows. In Sec. 2 we examine the breather 
lattice and two-soliton lattice solutions of the mKdV equation with positive sign of 
nonlinearity and discuss their stability. In Sec. 3 we examine periodic solutions of 
the mKdV equation with the negative sign of the nonlinearity. In Sec. 4 we briefly 
comment on the corresponding KdV solutions and follow that with our conclusions in 



Sec. 5. 



2. Breather and Soliton Lattices of the Positive mKdV Equation 



For a field u(x,i), the positive mKdV equation is given by 

u t + 6u 2 u x + u xxx = 0. 
Using the ansatz 

u = -2-^-[tan _1 0], 



(1) 



(2) 
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we are able to obtain several exact spatially periodic solutions. The first one is a 
breather lattice solution [T7| 

4> = asn(ax + bt + ao, k)dn(cx + dt + Co, m), (3) 

where ao,co are arbitrary constants. In fact, all the solutions discussed in this paper 
admit such constants even though we will not always display them hereafter. In this 
solution: 



c c 4 k b 

a' a 4 (1 — m) ' a 



[a 2 (l + k)-3c 2 (2-m)} 



- = [3a 2 (l + fc)-c 2 (2-m)], (4) 
c 

while sn(x, k), (cn(x, k) below), and dn(x, m) are Jacobi elliptic functions with modulus 
k and m, respectively. This solution was studied in detail previously in 17 . For 
completeness, we discuss briefly the relevant results. In order to ensure periodicity of 
the solution (in space), a commensurability condition was postulated in that work. In 
its strict form (strong commensurability) , this condition demands that the two elliptic 
function terms of Eq. (0) have the same period, i.e. AK{k)/a = 2K(m)/c (where 
K(k) denotes the complete elliptic integral of the hrst kind). In its weaker form (weak 
commensurability) , the condition demands that the periodicities arc rational multiples 
of each other i.e., 

ip m = 2q ^M (5) 

a c 

with p,q € Zi_, note that p = q = 1 yields the strong commensurability as a special 
case. Both strong and weak forms, however, resulted in unstable dynamical evolutions 
of the breather lattice |17) (induced by means of numerical perturbations to the exact 
solution). We further showed that this solution could be stabilized through ac driving 
and damping. In the limit k — > 0, m — > 1, but with k/(l — m) — v 4 , this solution goes 
over to the well known single breather (bion) solution 

<fi = —v sin(aa; + bt + ao)sech(cx + dt + Co). (6) 

Remarkably, it turns out that there is a different breather lattice solution to the 
positive mKdV equation given by 

4> = acn(ax + bt + ai, k)cn(cx + dt + ci, m), (7) 



where 



a 2 



mk c 4 fc(l — k) 

(1 — m)(l — k) a 4 m(l — m) 



- = [a 2 (l - 2k) + 3c 2 (l - 2m)}, - = [3a 2 (l - 2k) + c 2 (l - 2m)]. (8) 
a c 

Notice that in the limit k — > 0, m — * 1, a\ = ir/2 + gq c\ = cq but with k/(l — m) = v A , 
this solution also goes over to the same breather (bion) solution JHJ). However, for 
other values of k, m, the two breather lattice solutions © and J7J) are quite different 
(see also the comparison of the right panel of Fig. 0). We should note here that the 
characterization of the solutions that we present as breather or two-soliton lattices is 
given on the basis of what their limiting (see e.g., Eq. ©) profile looks like, as the 
elliptic functions asymptote to trigonometric/hyperbolic ones. 

The dynamical stability of the breather lattice solution 10 was examined for 
both strong and weak commensurability by means of direct numerical simulations. 
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Figure 1. The left panel shows the x-t (i.e., space-time) evolution of the contour 
of the unstable breather lattice solution of the form of Eqs. a = 1, k = 0.5, 

c = 1.707 and m = 0.03 (such that K(k)/a = 2K (m)/c) were used in the initial 
conditions. The middle panel shows the spatial profile at t = (initial condition) 
and at t = 55 (after the instability has set in). Finally, the right panel shows 
with a solid line the initial condition and with a dashed line the (other breather 
lattice) solution of Eqs. J3J-JIJ for the same parameters. 



The numerical scheme used here, motivated by the KdV discretization of Ref. [T§] 
was analyzed previously in |17| . However, the results were also verified with different 
discretizations of the integrable model such as the ones proposed in [20|- In Fig- Q] 
we demonstrate a typical example of the dynamical evolution of the solution of Eqs. 
Q-JHJ- We observe that similarly to the previously obtained breather lattice solution 
of the breather lattice identified above is dynamically unstable in the mKdV 
equation and results in a few more strongly and many weakly localized peaks. Hence, 
the instability of mKdV breather lattices appears to be generic. 

Using the ansatz in terms of Eq. J5J , we can obtain yet another novel family of 
solutions, namely a two-soliton lattice 

<j> = asc(ax + bt, k)dn(cx + dt, m), (9) 

where a = —c/a, sc(x, m) = sn(x, m)/cn(x, m) and 

4 = T— - = -[a 2 (2-fc)+3c 2 (2-m)], 
a 1 — m a 

- = -[c 2 (2-m) + 3a 2 (2-k)}. (10) 
c 

This solution can also be obtained from the breather lattice solution (J3J by taking 
a — > ia, b — * ib, a — * — ia and using the well known relations 

sn(x, 1 — m) dn(x, 1 — m) 

sn(ix, m) = l — - — , dn(ix, m) = — - — r- , 

cn(x, 1 — mj cn(x, 1 — mj 

cn(ix,m) = — - — \ -. (11) 

cn(x, 1 — m) 

In the limit k — > 1, m — * 1 but with (1 — k)/(l — m) = v A this solution goes over to 

<f> = -V sinh[a(x - [1 + 3u 2 ]i)]sech[av(x - [3 + v 2 ]t)] , (12) 

which, except for v = 1, is the well-known 2-soliton solution (hence the classification 
of this solution as a 2-soliton lattice). For v = 1, however, it corresponds to the 
one-soliton solution. The latter implies k — m and c = a from Eq. H10JI . 

An example of the dynamical evolution for the two-soliton lattice solution is shown 
in Fig. As it can be seen, this solution persists unchanged for long dynamical 
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Figure 2. The left panel shows the spatio-temporal evolution of the contour plot 
of a two-soliton lattice in the case of the positive mKdV equation for a = 0.25, 
k = 0.1, c = 1.596 and m = 0.999, such that K(k)/a = 2K(m)/c. The right panel 
shows the t = (top panel) and the t = 150 (bottom panel) spatial profiles. 



evolutions (i.e., for times of the order of 150 in the arbitrary time units of the time 
evolution of Fig. [SJ, hence our numerical simulations indicate that it is dynamically 
stable. 



3. Lattice Solutions of the Negative mKdV Equation 

The negative mKdV equation is given as 

lit - &u 2 u x + u xxx = . (13) 
In this case, to identify the corresponding solutions, we start with the ansatz 

u = -2— [taah _1 0l. (14) 
ax 

One can then show that the field 4> satisfies the equation 

(1 - (j> 2 )[(f>t + <t> xxx ] + &4>x[<t><t>xx - <t>l\ = 0. (15) 

Using an ansatz in terms of Jacobi elliptic functions we obtain the following new 
periodic solution: 

4> — adn(ax + bt, k)dn(cx + dt, m), (16) 

where 

4 = 7— - = -[« 2 (2-fc)+3c 2 (2-m)], 
a* 1 — to a 

o?= j= = . , - = -[c 2 (2-TO) + 3a 2 (2-fc)]. (17) 

Unfortunately, this solution may be singular for a finite x (for a given t) , depending on 
the parameters. In particular, the derivative of Eq. (|14f> induces a term ~ 1/(1 — u 2 ). 
However, from the properties of the elliptic functions, one can obtain that: 

A /(l-fc)(l-m)< M 2 <-== 1 (18) 

Vl 1 - fc )(! - m ) 
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Figure 3. Same as Fig. [5] but for the periodic solution of the negative 
mKdV equation with a = 1, fc = 0.5, c = 1.719 and m = 0.057 (such that 
K(k)/a = 2K(m)/c). The dynamical evolution is stable. 



Equation l|18ll . in turn, implies (given the continuity of u) that u will, typically, assume 
the value 1 for a certain x, hence that the corresponding solution will, gencrically, be 
singular. We thus do not consider it further here. 

Another, more interesting solution of the negative mKdV equation is given by 

4> = asn(ax + bt, k)sn(cx + dt, m), (19) 

where 

4 = -, - = [a 2 (l + £0 + 3c 2 (l + m)], 
m a 

a 2 = Vk^,, - = [c 2 (l + m) + 3a 2 (l + £;)]. (20) 

c 

For fc = to, this solution degenerates to the well-known soliton-lattice solution which 
in the limit fc = to = 1 goes over to the famous one-soliton solution of mKdV with 
negative sign. 

We have examined the solution (|19|l for strong, as well as weak commensurability. 
We have performed numerical simulations for various (p, q) pairs including the 
degenerate case with c = a amd fc = to. The principal feature of the temporal 
evolution of this solution lies in its dynamical stability, in the sense that for the 
duration of the numerical simulation (of the order of t = 150 time units) the structure 
robustly maintains its character, see Fig. 

4. Periodic Solution of the KdV Equation 

Through the Miura transform v = u 2 ± u x |18| we can obtain the periodic solutions 
of the KdV equation (v) from the corresponding solutions of the negative mKdV 
equation (u). In particular, to the solution (|19fl of negative mKdV, corresponds the 
KdV solution of the form: 



+ U *~ Z (1-0 2 )(1 ' ^ 



where prime denotes derivative with respect to the argument. The other solution 
u 2 — u'{x) is obtained from here simply by changing a to -a. Using the solution (|19|l 
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of the negative mKdV equation, we find that the solution jj^Jl can be expressed in the 
form 



Perhaps more interestingly, the Miura transform is an exact transformation 
between the solution of the negative mKdV equation (V(x,t)) and the one of the 
KdV equation. This implies that the numerically observed stability of the above 
mentioned lattice solution of the negative mKdV equation carries over to the existence 
and stability of such a solution in the setting of the KdV equation. 

5. Conclusion 

In this Letter, we have obtained new classes of periodic solutions for the positive 
and negative mKdV equations. The positive mKdV solutions J7J and © could 
be identified as the breather lattice and two-soliton lattice solutions, since in the 
appropriate limit they reduce to single breather and two-soliton solutions, respectively. 
In the case of the negative mKdV equation, lattice solutions have been identified in 
the form of Eqs. (|16|l and i|19fl . However, the latter have not been designated as 
breather or soliton lattices (as the process of obtaining limiting expressions is less 
straightforward for the negative mKdV case) . 

We have also examined the stability of the various periodic solutions of mKdV 
equation with both positive and negative signs. We have found that the two-soliton 
lattice solution of the positive mKdV and the lattice solution l|19|) of the negative 
mKdV are quite robust with respect to perturbations in contrast with the breather 
lattice solution of the positive mKdV equation. These results, and more specifically 
the stability of the negative mKdV lattice solution, have direct implications for the 
corresponding lattice solutions of the KdV equation. 

This work was supported in part by the U.S. Department of Energy. PGK is 
grateful to the Eppley Foundation for Research, the NSF-DMS-0204585 and the NSF- 
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